Disorder Effects in the Quantum Hall Effect of Graphene p-n Junctions 
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The quantum Hall effect in graphene p-n junctions is studied numerically with emphasis on the 
effect of disorder at the interface of two adjacent regions. Conductance plateaus are found to be 
attached to the intensity of the disorder, and are accompanied by universal conductance fluctuations 
in the bipolar regime, which is in good agreement with theoretical predictions of the random matrix 
theory on quantum chaotic cavities. The calculated Fano factors can be used in an experimental 
identification of the underlying transport character. 
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When a mono-layer of honeycomb lattice is sin- 
gled out of graphite,^ this two-dimensional mate- 
rial, dubbed graphene, acquires extraordinary elec- 
tronic properties j^i^i^ Electrons in graphene mimic mass- 
less Dirac fermions with extremely high mobility and 
tunability,^'^ which makes this material interesting both 
theoretically and practically. The tunability of the car- 
rier type via the electric- field effect, in particular, allows 
for the realization of graphene p-n junctions using only 
electrostatic gating. ^'^ The quantum Hall effect in these 
graphene p-n junctions has shown new fractional plateaus 
in the bipolar regime^ that were explained by uniform 
mixing among edge states at the junction interface.^ The 
mechanism of the mode mixing, however, is still unclear. 

In this paper we address this problem by investigat- 
ing the transport characteristics of graphene p-n junc- 
tions in the quantum Hall regime with disorder at the 
junction interface. Our calculations are based on the 
Landauer-Biittiker formalism for coherent transport, ^"^ 
and the results are explained using the random matrix 
theory (RMT) of quantum transport F^iii^ 



region 1 




FIG. 1: (Color online) Schematic diagram of a graphene junc- 
tion. Two locally gate-controlled regions of a graphene strip 
(width W), connected with a reservoir at each of the far ends, 
are jointed by a transition area (length d) where potential is 
assumed to be composed of a slope along the strip and a ran- 
dom distribution over each site within. The lower curve shows 



the profile of the relative Fermi energy E' 



(r) 



The setup of our simulation is illustrated in Fig. [T] 
A graphene strip of width W is divided into two regions 
by a transition area with length d. In either region the 
carrier type and density can be locally tuned through an 
electrostatic gating, and the relative value of the Fermi 
energy to the local charge neutrality point is defined as 

(r) 

the relative Fermi energy E\ ' . The whole sample is sub- 
ject to a perpendicular magnetic field, and the Landau 
levels are formed in the quantum Hall regimci^ With a 
specific Fermi energy, the filling factor of Landau levels 
is scarcely dependent on the details of the sample edge, 
provided the sample size is big enough^i^ therefore we 
use samples with zigzag edges to carry out our simula- 
tion but the results are applicable to general cases. On 
the energy scale of our problem, both Zeeman splitting 
and spin-orbit interaction are negligibly small, different 
spin states can be taken as degenerate and uncorrelated, 
thus we assume spin is irrelevant in our calculation and 
simply multiply the result by a factor accounting for the 
spin degree of freedom. For this reason the filling factors 
vi and 1^2 given below are all "spinless", i.e., their values 
are ±1, ±3, ±5 • • • instead of ±2, ±6, ±10 • • • . 

Magneto-conductance of a graphene p-n junction has 
been theoretically discussed in terms of its valley isospin 
dependence.— These discussions in general assume an ab- 
sence of intervalley scattering across the graphene p-n 
junction. To explain the fractionally quantized plateaus 
observed in the experiments,^'^ however, full mode mix- 
ing was required in random matrix theory, which in- 
evitably involves intervalley scattering^. We attribute 
the inter-valley scattering in our model to disorder po- 
tential which varies quickly enough in the scale of the 
lattice constant. Disorder in graphene may have various 
sources, and the understanding of its role in transport 
properties is still incomplete^ii^ 

In this paper we focus on the effect of the random po- 
tential at the interface of the junction. The reason is 
that current inside either region 1 or region 2 is carried 
by quantum Hall edge stateSfii thus immune to most dis- 



order inside either region, while scattering among states 
at the interface of the two regions contributes to the ma- 
jor effect of disorder on the overaU conductance. Disor- 
der at the interface may come from intrinsic sources hke 
vacancies and impurities, or from extrinsic sources hke 
random potential introduced by the irregularities of the 
gate edge. We model the disorder potential at the inter- 
face by using the Anderson-type on-site energyiS which 
varies randomly from site to site within the transition 
area, where a potential slope connecting two sides of the 
junction serves as the background potential. Hence the 
total Hamiltonian reads 



i/ = Ve(r, 
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where c] and Ci are the electron creation and annihilation 
operators at site r^ = {xi,yi), respectively, t w 2.8 eV 
is the nearest neighbor hopping energy in the graphene 
lattice and {i,j) stands for a nearest-neighboring pair. 
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is the phase acquired when an electron hopping from rj 
to Ti in an external field B described by vector potential 
A, and the on-site energy £{ri) are ei and £2 m region 1 
and 2, respectively; 

X- 1 
£(n) = {£2 - £i)-j + 2 (ei + £2) + rand(A) (3) 

in the transition area {\xi\ < d/2) with rand(A) a ran- 
dom number uniformly distributed in [—A/2, A/2]. The 
width of the junction in our simulation is taken to be 
W = 200(V3a/2) « 42.6 nm, where a « 0.246 nm is 
the lattice constant of graphene, and the length of the 
interface area is taken to be d = 20a « 4.9 nm. Landau 
gauge A{r) = {—By, 0) is adopted, with the magnitude 
of the magnetic field S w 113 Tesla which is equivalent 
to a magnetic flux $ = h/701e in each unit cell. The 
magnetic length in this case is Ib — yWi^B) « 2.4 nm, 
which is about a half of the length of the interface area, 
or 1/18 of the width of the sample. It should be men- 
tioned here that in a real sample which is presumably 
much larger in size, the magnetic field necessary for the 
quantum Hall effect can be much smaller. 

We assume coherent transport in the graphene junc- 
tion, where the Landauer-Biittiker formalism can be 
applied>i2i Transmission functions Tpq {p,q — 1,2 and 
p ^ q) oi the junction are calculated by using the recur- 
sive Green's function technique.— i^ Vanishing net cur- 
rent in equilibrium implies that T21 = T12 = T, therefore 
the conductance is proportional to either of the transmis- 
sion functions, G = (e^//i)T, where the spin degeneracy 
has been included in T, and the variance of the con- 
ductance Var(G') = (e2//i)2Var(T), where Var(T) rep- 
resents the variance of T. In the following we will be 
satisfied with observing only the behavior of T in dif- 
ferent situations. Each situation, that is, each experi- 
mental condition under which measurements are made, 
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FIG. 2: (Color online) (a) The mean and (b) the variance of 
the transmission function T as functions of the relative Fermi 
energy Ef in region 2, with fixed relative Fermi energy E^ 
in region 1 and fixed disorder strength A. Cases with difiterent 



pirl) 



(rl) 



Ef and different A are also compared. E, > implies 



r('-i) 



n-type of region 1, Er jt — 0.10 corresponds to the filling 



factor vi — 1 and Ef jt — 0.15 corresponds to v\ 



is identified with a specified combination of £1, £2 and 
A, while various configurations of disorder are subject to 
some self-averaging process in each measurement. This 
self-averaging process could be a result of time dependent 
electric field used in the experiments,— and will suppress 
the fluctuation of the measured conductance, thus makes 
the mean value a reasonable account for the experimen- 
tal observation. Our calculation extracts the mean and 
the variance of the transmission functions T in each situa- 
tion from output of 40,000 samples with different disorder 
configurations. 



III. RESULTS 

The calculated transmission functions as shown in Fig. 
[2t^a) have surprisingly recovered the quantized transport 
plateaus observed in the experiment by Williams et alii 
In junctions with the disorder strength A = 2< (solid 
lines in Fig. [5]), the ensemble average of the transmis- 
sion functions form nearly perfect plateaus in the bipolar 



regime (FTt < 0), and the height of each plateau is 



(r) = 2 
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(4) 



with j/totai = \v\\ + \v2\- Corresponding to each plateau 
of the averaged transmission function, the ensemble vari- 
ance of T also develops into a plateau described by 



Var(T) = 4 X 



{V\V2Y 



(kii + i^2ir[(i^ii + 1^21)^-1] 



(5) 



1.5 



0.0 

0.4 

0.3 

ro 0.2 
> 

0.1 
0.0 



1 




^ 






*^> 


\ 


v,=-l 

V2=-3 


(a) 


^1 


L 








(b)"" 


..^ 


/\ 






/ 


x 






^s^ 



5/3 
3/2 



3/20 
5/63 



4 6 

A/t 



10 



FIG. 3: (Color online) (a) The mean and (b) the variance 
of the transmission function T as functions of the disorder 
strength A in bipolar junctions. Region 1 is of n-type with 
the filling factor z^i = 1, and region 2 is of p-type with the 
filling factor V2 = —1, —3, —5, respectively. The shadow high- 
lights the regime where both (T) and Var(r) show plateaus 
predicted by Eq. (|4]) and Eq. ([5]), respectively. 



Both Eq. (ID) and Eq. (O are the predictions of the 
RMT on a quantum chaotic cavityfi2, with the additional 
factors 2 and 4 from the spin degeneracy. In the unipolar 
regime [Et > 0), plateaus of the ensemble average (T) 
are only partly formed when A = 2t, and the height may 
not be accurately of the expected values given by 



(T) = 2xmin(|z.i|,|i/2|) 



(6) 



The transmission functions show large ensemble variance 
where (T) has a large deviation from Eq. ©. Decreased 
disorder strength in the junction interface, in contrast, 
leads to better-developed plateaus of (T) in the unipolar 
regime, at the cost of losing the quantized values of (T) 
and Var(r) in the bipolar regime. This is presented as 
the dash lines in Fig. [5] with A = 0.5t. 

The plateaus described by Eq. (|4|) in the bipolar 
regime and described by Eq. ^ in the unipolar regime 
are the signature of the quantum Hall effect in a single 
graphene junction,^ though in some cases the accuracy of 
the plateau is poor in the experimental dataf^i^ compared 
with the expected value. By taking into account disor- 
der in the junction interface area and a self-averaging 
process in the measurement, our calculations clearly pro- 
duce these plateaus. In addition, the lack of the accuracy 
of the plateau height is attached to the strength of the 
disorder, and is reflected in the variance of the transmis- 
sion functions. 

The experimentally observed conductance plateaus of 
a bipolar graphene junction in the quantum Hall regime 
have been explained as the result of the complete mix- 
ing of quantum Hall edge states at the junction inter- 



face due to scattering, and the departures of the exper- 
imental data from Eq. Q have been attributed to the 
incomplete mixing of edge statesiii^ We emphasize here 
that because the spin-flip process is negligible in this sys- 
tem, the mode mixing can only happen among states of 
the same spin quantum number. Thus to correctly ex- 
press (T) and Var(T) in terms of the filling factors i^i 
and 1/2, the filling factors must be spinless, with the spin 
degree of freedom included as an extra multiplier to T. 
Such expressions, compared with the expressions using 
the spinful filling factors ±2, ±6, ±10...)^ show no quan- 
titative difference as for (T), but significant differences 
as for Var(T). 

The disorder dependence of the transmission functions 
with specific combinations of filling factors in the bipolar 
regime is shown in Fig. [31 The shadowed region (roughly 
1.7t < A < 2.8i) is where complete mode mixing hap- 
pens in all three cases, i.e. i^i = 1 and 1/2 = —1, —3, —5, 
respectively. The averaged transmission functions de- 
velop into plateaus of height described by Eq. (|4]) si- 
multaneously in this region, and the ensemble variances 
of the transmission functions also develop into plateaus 
predicted by Eq. dS]). In the language of the RMT^ii^ 
the ensembles of scattering matrices S (of a specific spin 
quantum number) under these circumstances are the cir- 
cular unitary ensembles (CUEs), that is, S matrices in 
these ensembles are uniformly distributed over the uni- 
tary group Z/^(t'totai)- Average over the CUE is equal to 
an integration over the unitary group. And it is this inte- 
gration that lead to the " universal" value of the averaged 
transmission function Eq. ^ , and the universal conduc- 
tance fluctuation (UCF) given by Eq. ([5]). The graphene 
bipolar junctions in this disorder regime (shadowed) are 
nearly ideal realizations of the quantum chaotic cavities 
characterized by Eq. ^ and ([5]) . 

The ensemble of S matrices, however, is actually de- 
pendent on the disorder strength A, which represents 
how much the scattering potential can be varied from 
one sample to another. In Fig. [3] we see that to the left 
of the shadowed region (roughly A < l.7t), both (T) and 
Var(T) deviate from Eq. ^ and ([5]) with decreased A, 
indicating the deviation of the actual ensembles of 5* ma- 
trices from the CUEs. In other words, these are the cases 
of incomplete mode mixing. Notably, when i/i = — 1/2 = 1 
there is an extended range of A where the plateau of (T) 
is preserved. This is easily understood as there are only 
two modes to be mixed in this case. When the disor- 
der at the junction interface is larger than the range of 
the shadowed region in Fig. [31 on the other hand, the 
conducting edge modes in either region 1 or region 2 are 
expelled from the junction interface area gradually, scat- 
tering between modes from one side of the junction to 
the other is weak, both (T) and Var(T) decrease until 
the whole junction is "cut off' . This is the case of what 
we see in the large-A regime of Fig. [3] 

Compared with the bipolar junctions, transport in a 
unipolar junction with disordered junction interface has 
a more straightforward picture. The conductance of a 




FIG. 4: (Color online) The averaged transmission function 
(T) as a function of the disorder strength A in unipolar (n-n) 
junctions with different combinations of filling factors. The 
inset shows the corresponding variance of T. The regime is 
shadowed where all plateaus of (T) are preserved. 
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FIG. 5: (Color online) The Fano factors in the bipolar regime, 
as functions of the relative Fermi energy E)!' ', with ui = 
f (black) and 3 (red/dark gray), respectively, and A — 2t. 
Plateaus predicted by Eq. ((Tjl are indicated by short dash 
lines. 



Fano factors develop into plateaus described by 



unipolar junction is mainly contributed from the compat- 
ible edge modes in the two regions, and the tunneling of 
carriers from one region to another mainly happens near 
the lateral edges. The disorder at the interface impedes 
the coupling of these modes to reduce the transmission 
through the junction, which is contrary to the bipolar 
cases. The calculated (T) and Var(r) in the unipolar 
regime are shown in Fig. [4] with different combinations 
of filling factors. In each case (T) exhibits a quantized 
plateau when the disorder is weak (shadowed), and starts 
to deviate from the plateau at a critical value of the dis- 
order strength where Var(T) also begins to deviate from 
zero. We notice that the shadowed regions in Fig. [3]and[3] 
do not overlap over the range of A, which implies that (T) 
do not develop into plateaus simultaneously in the bipo- 
lar regime and the unipolar regime, especially when the 
filling factors are large. This fact highlights the opposite 
roles the disorder plays in the formation of conductance 
plateaus in the bipolar and the unipolar junctions. It sta- 
bilizes the plateaus in the bipolar case when its intensity 
is in a certain nonzero range, while it tends to destroy the 
plateaus in the unipolar case at the same time, though 
in the limit of strong disorder electrons will be blocked 
from tunneling through the junction in both cases. 

The character of the quantized transport in these junc- 
tions can be experimentally identified by measuring the 
electron shot noise. ^^ The Fano factor, defined as the 
ratio of the actual shot noise to the Poisson noise, is 
extracted from our simulation using the equation F — 
(^ Tn{l — Tn)) / (Y^ Tn), whcrc the summations are taken 
over different transmission eigenvalues indexed by n. It 
is found that the Fano factors corresponding to the con- 
ductance plateaus in the unipolar regime are identically 
zero, as expected from the dissipationless transport via 
the quantum Hall edge modes. In the bipolar regime, the 
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(7) 



corresponding to the plateaus of (T) , as shown in Fig. [5l 
Eq. ([7]) is again a straightforward outcome of the RMT 
applied to a quantum chaotic cavity with few transmis- 
sion modesi^ It is quantitatively different from the Fano 
factors discussed in reference^, especially when the num- 
ber of the transmission modes is small. And this is a 
key to examine experimentally the existence of the UCF 
revealed in our simulation. 



IV. DISCUSSIONS AND CONCLUSION 

Before ending this paper we address the issues on the 
sources and types of disorder in graphene which is rel- 
evant to this study. The experimentally observed con- 
ductance plateaus were explained as a result of the full 
mixing of these quantum Hall edge modesi^i^ The details 
of the mode mixing are determined by the form of the dis- 
order and in turn provide information about the disorder 
therein. The quantum Hall edge modes in a graphene 
strip can be indexed by two quantum numbers, spin and 
valley isospiuj^ besides the Landau levels they belong to. 
Thus the mechanism of the mode mixing inevitably in- 
volves the presence of spin-flip scattering and/or interval- 
ley scattering. Considering the negligible magnetic impu- 
rities and spin-related interaction in the current graphene 
samples, and ignoring the magnetic edge states reported 
in some specific graphene ribbonsjS4, we assume the ab- 
sence of spin-flip scattering in this study, and use the 
"spinless" filling factors i^i and 1^2 (= ±1, ±3, ±5 • • • ) in 
our expressions 1^ The presence of intervalley scattering 
among the quantum Hall edge modes in the graphene p- 
n junctions, however, is still a question. Previous works 



on the disorder effects in graphene p-n junctions gener- 
ally assumed weak intervalley scattering, and discussed 
the valley-isospin dependence of the conductance.— For 
example, in the work by Tworzydlo et al.^ it was as- 
sumed that each impurity has the Gaussian potential 
profile C/iexp(— |r — iiip/2^^) of range ^ and random 
height Ui G {—5,6). A large ^ implies a long range of 
the disorder potential which suppresses the intervalley 
scattering, while a small £, implies (5-function-like disor- 
der potential which provokes the intervalley scattering. 
The former case was studied but the plateaus observed 
in the experiments were not properly recovered. That is 
the reason why we choose the latter case as our starting 
point in this work, though the origin of the short-range 
disorder potential is not completely understood. It turns 
out our calculations do produce the reported conduc- 
tance plateaus as in Eq. ([4]) , as well as quantities like the 
Fano factors as in Eq. ([7]) that can be further examined 
experimentally. Still it is interesting to investigate the 
crossover effects between long-range disorder potential 
and short-range disorder potential in the junctions, and 
the results would possibly explain the not-fully-developed 
conductance plateaus in the experiments /'*■ This investi- 
gation will be a subject of our later work. Furthermore, 
there also exist other forms of disorder that cannot be 
described by random local potentialj^ii^. For example, 
effective gauge fields can be induced by ripples, topo- 
logical lattice defects, strains, or curvatures etc. in the 
graphene sample, and will affect the Landau levels, es- 
pecially the lowest one^ We believe that these various 
forms of disorder will affect the behavior of the edge- 
mode mixing in various manners, but this is out of the 
scope of the current work. 

In short, our numerical simulation of the quantum Hall 
effect in graphene p-n junctions has reproduced the quan- 
tized conductance plateaus observed in the experiment. 
The UCF and quantized values of the Fano factors are 
found to be accompanying the conductance plateaus in 
the bipolar regime, which is well explained by the RMT 
of quantum transport. The bipolar graphene junction 
in the quantum Hall regime mimics an ideal quantum 
chaotic cavity, which is another example of the extraor- 
dinary transport character of graphene. 
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V. APPENDIX: DERIVATION OF EQS. (gD, ©, 
AND dlD 

In this appendix, we derive Eqs. dl]), ©, and ^ using 
the random matrix theory. 



Random matrix theory of quantum transport is based 
on the assumption that the scattering matrix S" of a 
chaotic cavity is uniformly distributed over a specific 
group (the so-called circular ensemble) determined by the 
symmetry of the Hamiltonian. Basically there are three 
classes of the groups.^* If time-reversal symmetry is bro- 
ken ((3 = 2), S is only constrained by unitarity, which is a 
result of current conservation, thus S belongs to a unitary 
group. If time-reversal symmetry is preserved together 
with the presence of spin-rotation symmetry (/3 = 1), 
then 5' is both unitary and symmetric: S = S'^ , where 
the superscript T indicates the transpose of the matrix, 
this leads to an orthogonal group. If time-reversal sym- 
metry is preserved but spin-rotation symmetry is broken 
(/3 = 4), which is the case when spin-orbit interaction is 
present, then S is unitary and self-dual: S = S^ , where 
the superscript R indicates the dual of a quaternion ma- 
trix, the group is called symplectic. 

In random matrix theory, the statistics of transport 
properties is obtained from the statistics of an appro- 
priate circular ensemble. For example, the mean of the 
transmission probability T„„i = | S',„„ | ^ is given by 



{Tn 



dll[b)Snmi>nm 



(8) 



where n and m stand for transmission eigen-channels. 
We will omit the mathematical details of formulating the 
measure dfj,{S) in a group space that S belongs to, and 
use the following two equations as being established)^ 



{Ua.aU, 
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T^afsSab 



(9) 



{UaaUa'a'U^fjUp/f,,) CUE 
= J;H Ti^afjSabSa'P'Sa'b' + Sap' Sab' Sa' pSa'b) 

~~nTn^ 7T('5Q/3^ab'^a'/3"5a'6 + Sap' SabSa' f3Sa'b') (10) 

where U is an NxN unitary matrix belonging to the 
circular unitary ensemble (CUE). Since the system we 
are going to discuss only involves the CUE, we will drop 
this notation hereafter. 

Starting from these equations, the conductance (G), 
its fiuctuation Var(G), and the Fano factor F for the 
fully-mixed quantum hall edge transport in a graphene 
p-n junction can be readily calculated. We first assume 
no spin-flip scattering is allowed, so that S matrix for 
each spin component is unitary by itself, the case includ- 
ing spin-flip scattering will be discussed in what follows. 
Suppose there are 2Ni and 2N2 edge states (they are also 
eigen-channels away from the junction area) at the two 
sides of the junction, respectively, the full S matrix is of 
dimension 2{Ni + N2), but it is divided into two uncou- 
pled spin subspaces in the absence of spin-flip scatter- 
ing, therefore matrices in the CUE are only of dimension 
A^totai = A^i -t-iV2 in terms of each spin subspace. Also the 



strict spin degeneracy implies that S'^„^ — S',„^ = Snm- 
Regarding this we have 



e 



Ni N2 

Ni N2 2 

n=l m—1 

2N1N2 e^ 
Ni + N2~h' 



(11) 



In contrast to the spin-conserved scattering case, the 
mixing among different spin modes will lead to quantita- 
tively different variance and Fane factors, while leaving 
the conductance mean unchanged. In the random ma- 
trix theory, this is derived as following. Still suppose 
there are 2Ni and 2N2 edge states at the two sides of the 
junction respectively, the full S matrix is still of dimen- 
sion 2(A^i -I- A^2), but the spin-off-diagonal parts are now 
non-zero, and the 5' matrix is unitary only as a whole. It 
is straightforward to show that in this case the formulae 
corresponding of Eqs. ([1]), (O, and ([7]) are given by 



Var(G) = {(G^) ~ {Gf 

I / Ni N2 

~ \ \ / J / J / y -'-nrn 
\ \{T— ±1 n— 1 m—1 
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n,n'=^l m,?n' — 1 



4N^N^ 



{Ni + N2y [{Ni + N2)^ - 1] 



-{G? 



(12) 



The Fano factor F = (Trace(ttt(i _ ttt)))/(Trace(tt^)), 
where t is the transmission submatrix of S,^^ and 
Trace(ttt) is nothing but G/^. Thus 

iV]__KA^ >r^ V^ ^ iq<y qa*o<y o<t * \ 



2N1N2 
N1N2 



cr— ±1 n,n' — 1 m.ni' — l 



{Ni + N2Y - 1 



(13) 



In this way we have derived the formulae expressed in 
Eqs. dH), (O, and ([7]) in the framework of the random 
matrix theory. 



(G) 



Var(G) 



and 



4^1 iV2 e2 
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{2Ni + 2N2Y [{2Ni + 2N2Y ~ 

[Ni + N2f [{Ni + N2Y - 1/4] U ' ^ ^ 



F 



AN1N2 



{2Ni + 2N2)- - 1 
(7Vi+iV2)'-l/4' 



(16) 



respectively. Except for a coincident equivalence of the 
mean values of the conductance (G), the quantitative 
differences of the UCF Var(G) and the Fano factors F 
between two cases, namely whether scattering is spin- 
conserved or not, are obvious, especially when the values 
of A^i and A''2 are small. 
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